Phase ordering of two-dimensional symmetric binary fluids: a droplet scaling state 
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The late-stage phase ordering, in d = 2 dimensions, of sym- 
metric fluid mixtures violates dynamical scaling. We show 
however that, even at 50/50 volume fractions, if an asym- 
metric droplet morphology is initially present then this sus- 
tains itself, throughout the viscous hydrodynamic regime, by 
a 'coalescence-induced coalescence' mechanism. Scaling is re- 
covered (with length scale I ~ t, as in d = 3). The crossover 
to the inertial hydrodynamic regime is delayed even longer 
than in d — 3; on entering it, full symmetry is finally restored 
and we find / ~ t 2 ^ 3 , regardless of the initial state. 

When a symmetric binary fluid is quenched below its 
critical temperature it will phase separate into domains. 
This initial stage of phase separation is diffusive; but in 
this letter we are concerned with the hydrodynamic stage 
of the process, which follows afterwards. The evolution 
equations can be taken as the Navier- Stokes and the 
advective diffusion equations: 

pD t u = - VP - (jNp^ + r?V 2 u (1) 
D t </> = V (MV W ) (2) 

Here u is the fluid velocity, D t the material derivative, P 
pressure, the chemical potential for the order param- 
eter 0, r] the viscosity, p the density and M an Onsager 
coefficient. As usual, \i§ = 5F/5<j> with F[<j>] a (symmet- 
ric) free energy functional; once the interfaces are sharp, 
the source term (~ </>V/i</>) in Eq. [l] is just the Laplace 
pressure of interfacial curvature |l[ . 

These equations are nonlinear, but deterministic; ther- 
mal noise should be irrelevant at late times. The only 
source of noise is then in the initial condition, and it is of- 
ten assumed that the physics is the same for all conditions 
without long-range correlations . But we find below 
that the nonlinear dynamics are more subtle: for two 
dimensional 50/50 binary fluid mixtures of equal viscos- 
ity, the hydrodynamic coarsening depends very strongly 
on the topology [Q of the sharp interfaces formed by 
the early time (short length-scale) physics, previously 
assumed irrelevant. Specifically, droplet morphologies 
evolve quite unlike bicontinuous ones: dynamical scal- 
ing, absent for bicontinuous states in d = 2 dimensions 
H, is recovered for a new, self-sustaining droplet state. 

Such scaling requires that any averaged property of 
the system is invariant if all lengths are scaled by a single 
length scale L(t). Whereas at early times diffusive phase- 
ordering leads to L oc t 1 ! 3 in both d = 2 and d = 3, scal- 
ing of interfacial curvatures with L implies that at later 
times the chemical potential obeys p^ ~ <j/L(t) where a 
is the surface tension (computable from F [</>]). Dimen- 
sional analysis of Eqs.1,2 then identifies two subsequent 



regimes, the viscous hydrodynamic (VH |6|]) and the in- 
ertial hydrodynamic (IH [Q) in which the driving force 
is balanced by viscosity and inertia respectively [|| . 

There is good numerical evidence for this behav- 
ior in bicontinuous d = 3 mixtures Q: allowing for 
a nonuniversal early time diffusive offset, and choos- 
ing the reduced unit of time as to = rf / pa 2 = 1, 
one finds a universal scaling curve l(t) where l(t) = 
L(t)pa/r] 2 is the reduced coarsening length. Choosing 
L = 2tt J S(k)dk/ J kS(k)dk with S(k) the density au- 
tocorrelator, one then finds that the viscous and inertial 
asymptotes are I = 0.07t and I = 0.9t 2 / 3 ; these cross at 
t* ~ 10 4 , a large value (though formally still 0(1)) §]. 

In two dimensions, in contrast, the VH regime was 
recently found to violate scaling altogether [|]]|]. Start- 
ing from a 50/50 bicontinuous state, a cascade of small 
droplets is left behind by the large connected domains as 
they coarsen; these droplets, which must coalesce in or- 
der to grow, evolve more slowly than the larger domains 
around them. Length measures defined (e.g.) from dif- 
ferent moments of the structure factor do not scale. Fu- 
rukawa ^ found, however, that the size of the largest 
domains does scale as I ~ t throughout VH; and that, in 
the IH regime, full scaling is recovered with I ~ t 2 ! 3 . 

For symmetric binary fluids, droplet morphologies usu- 
ally arise for an asymmetric composition. Until recently 
these morphologies were held to give non-hydrodynamic 
phase ordering by either the evaporation-condensation 
mechanism or the Brownian coalescence of droplets, each 
of which has t 1 ! 3 growth [Q. But Tanaka Q discovered 
that such droplet states can coarsen via a 'coalescence 
induced coalescence' (CIC) mechanism; Nicolayev et al. 
[ p"H considered this for some simplified d = 3 geometries, 
and predicted the recovery of a VH (I ~ t) growth regime 
mediated by CIC for </> > 0.26 in d = 3. 

In d = 2, droplet morphologies can, however, also 
arise in symmetric (50/50) quenches whenever there is 
asymmetry in the dynamics. This occurs, e.g., when the 
two components have different diffusivities so that M is 
strongly ^-dependent ]l2| ; or when one phase is viscoelas- 
tic [^3|. In the first case, the evolution equation should 
approach that of a fully symmetric mixture on entering 
VH. In the second case (since the mean shear rate is of or- 
der Ijl ~ 1/t), one must wait until t exceeds the Maxwell 
time; this could still be early in VH [Q. In either case 
the dynamic asymmetry is then forgotten, except that it 
has instilled local correlations in the "initial" condition 
for subsequent events. It is usual to argue that these 
correlations cannot influence the late time scaling states 
once the domain scale L(t) is larger than the original cor- 
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relation length; indeed, Bhattacharya et al. || thereby 
argued that viscoelastic phase separation could not lead 
to new scaling states. We now show that, in contrast to 
such expectations, a distinct CIC droplet scaling state 
does exist for a 50/50 quench, and survives throughout 
the VH regime. 

We use a well-established lattice Boltzmann method 
[|l5f to simulate Eqs. 1,2. We initialize each simulation 
with a droplet morphology, at 50/50 volume fraction, 
generated from a white-noise initial state using a separate 
viscoelastic lattice Boltzmann model [[l3|. (The initial 
droplets are at rest.) Thereafter the evolution equations 
are fully symmetric: the two fluids are Newtonian and 
have equal mobilities and viscosities. We choose param- 
eters such that the phase ordering is dominantly hydro- 
dynamic, ensuring that residual diffusion is negligible at 
the < 1% level Q|. The algorithm shows good isotropy 
for the shape of an equilibrium droplet ]l6| , but some dy- 
namic anisotropy is noticeable in the IH regime (visible 
in Figure 2). Reduced physical units l,t are defined as 
above, and nonuniversal (diffusive) time offsets for each 
simulation run calculated Q. All the simulations shown 
here were performed on a 1024 2 lattice, but terminated 
at L < 256 to mininize finite size artifacts The sur- 
face tension is 10 -2 , the density is 2, and we varied the 
viscosity from 0.5 to 0.001 (all in lattice units). 

We used three different length-scale measures [§,0, 
derived from the number of domains TV, the interface 
length Li and the derivatives of the order-parameter </>: 
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where A is the average of the eigenvalues of the derivative 
matrix <9 x </><9 x </>, and y are the sample dimensions. The 
lengths and are normalized such that each 

gives the unit cell spacing for a square lattice of droplets 
of 50% volume fraction. For bicontinuous morphologies 
in d = 3, all are comparable to L(t) as defined previously; 
for example, L d c± 1.2L. 

We show in Fig. 1 the evolution of the droplet mor- 
phology for a simulation run that lies wholly within the 
VH regime. After a brief transient (leading to the first 
coalescence) the coarsening mechanism is CIC, as first 
seen in experiments by Tanaka |l(| . The results are rem- 
iniscent of the predictions for CIC in d = 3 ]Tl| ], where 
it was predicted that the droplets that just underwent 
a coalescence would be the next to coalesce. Instead we 
observe that in d = 2 the coalescence of two droplets usu- 
ally induces a coalescence elsewhere. (This is especially 
clear in animations [p8|.) Perhaps for this reason, the 
CIC process does not (in d = 2) lead rapidly to percola- 
tion of the droplets as was suggested for d = 3 |11|. Such 
percolation would cause a rapid crossover to bicontinu- 
ity. Instead, the system achieves a self-sustaining droplet 
state that, in contrast to the bicontinuous analogue [[|, 



shows good dynamical scaling (Fig.l). This behavior is 
seen for all runs that lie within the VH regime as identi- 
fied from the universal scaling plot found below (Fig. 3). 
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FIG. 1. (a) The droplet scaling state for r\ = 0.22 (VH 
regime) for reduced times t = 1100, 1500 and 2400 respec- 
tively, (b) A section scaled by l(t) is shown; note the scale 
invariance. (c) The evolution of L # , L d and L° with time 
(where L # < L d < L°). All three length-scales agree well 
with each other and approach L ~ t (grey line). 

This behavior is quite unexpected [|| and only partly 
explicable. On one hand, while only local rearrangements 
(of order parameter or interface) are needed to transform 
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one to the other, droplet and bicontinuous morpholo- 
gies are very different: the largest connected domains 
present are small in one case, and percolate in the other. 
The issue is whether the 'chain-reaction' of CIC is fast 
enough to percolate the droplets [[Tl| : we find it is not fast 
enough. On the other hand, it is just fast enough to keep 
going indefinitely, whereas it might equally have stopped 
in finite time, leaving classical, non-hydrodynamic coars- 
ening processes to take over. Note that under the CIC 
process, large domains can no longer coarsen simply by 
deforming their tortuous boundaries, but have to coa- 
lesce. This explains why scaling is recovered: all domains 
are droplets and coarsen via the same mechanism. 

In Fig. 2 we see the counterpart of Fig.l for a run that 
lies wholly in IH. Here the initial droplet morphology 
(not shown) crosses over rapidly to a bicontinuous scal- 
ing state seen previously for the IH regime in two dimen- 
sions [§,[)). ^he mechanism for the rapid morphological 
transition is fluid inertia: two coalescing droplets wobble 
violently as they fuse together, causing local as well as 
distant coalescence events and rapid percolation of the 
initial droplet structure. In the IH scaling state, the dif- 
ference between droplet or bicontinuous initial conditions 
has thus been lost and local correlations inherited from 
early-time asymmetry do not have any long-term effect. 
The observed coarsening mechanism, in which the local 
topology is constantly being ruptured and reconnected by 
under damped inter facial oscillations, is somewhat remi- 
niscent of the 'turbulent remixing' process hypothesized 
in [|l9[ (see also [po| ). (This is somewhat hard to ascer- 
tain from the still pictures in Figure 2 but is easily seen 
in our animations fig] .) But in contrast to the prediction 
|L9| that I < CHt 1 * 7 ), we see (as in g) a simple I ~ t 2 ^ 
scaling throughout the accessible IH regime. 

Just as in d = 3, we cannot rule out a further crossover 
to a regime where the viscosity is so low that the growth 
rate is not limited by the input of energy from the surface 
tension, so that the primary balance in Eos. |l],|| is be- 
tween inert ial and viscous effects. Kendon || considered 
this case in d = 3 and argued that I ~ t 2 / 3 might survive 
for the structural length scale, so long as other lengths 
(describing the velocity field) do not scale the same way. 
It would be interesting to pursue this suggestion in d = 2, 
although the properties of d = 2 and d = 3 turbulence 
pl[ are very different. 

To establish scaling with confidence, Figure 1(c) and 
2(c) are not sufficient; we must combine several simula- 
tions with different viscosities to span many decades on 
the l(t) curve As shown in Fig. 3, this yields very 
good data collapse onto a universal scaling curve for all 
three length measures. A prefactor deviation appears be- 
tween and the others in the IH regime, but this results 
from the change in morphology (Fig. 2), not a violation 
of scaling. The asymptotes are (for lengths defined via 
L° and L d ) I = 0.02t and I = 0.9t 2 / 3 . These cross at 
t = t* ~ 2 x 10 5 , a value even larger than the d = 3 case, 
whose asymptotes (using L = 2tt J S(k)dk/ J kS(k)dk) 
are also shown for comparison. 
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FIG. 2. Recovery of bicontinuous scaling state for 
r\ — 0.0044 (IH regime), (a) The symmetric scaling state 
that rapidly evolved from the initial droplet morphology for 
t = 2x 10 7 , 4.2 x 10 7 and 9.3 x 10 7 . (b) Images scaled by /(£), 
showing dynamic scaling, (c) L° ,L d and L # all scale as t 2 ^ 3 
but L* is now larger that the other two measures, indicating 
the change in morphology. 

In comparable units, our VH droplet growth is signifi- 
cantly slower than for the d = 3 bicontinuous VH scaling 
morphology. This is reasonable: the amount of interface 
actually driving the coarsening is much smaller since only 
droplets undergoing coalescence contribute to the flow at 
any time, while in the bicontinuous case the entire in- 
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terface is in motion. Once the IH regime is entered the 
d = 2, 3 structures are both bicontinuous and the coeffi- 
cients of t 2 / 3 become much more similar. 
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FIG. 3. Universal scaling curve l(t) for the evolving 
droplet state. Note the separation of L # (on top) from L° and 
L d at the crossover to the IH regime, indicating the change 
in morphology. Asymptotes in grey; those for bicontinuous 
d = 3 coarsening Q are also shown. 

To summarize, in this letter we have found an unex- 
pected droplet scaling state in the phase-ordering of two- 
dimensional symmetric binary mixtures at 50/50 com- 
positions. This implies that the local order parameter 
correlations brought about by asymmetric early-time dy- 
namics (e.g. diffusion or viscoelasticity) can continue to 
control the behavior even when the coarsening length has 
become very large compared to the original length scale 
of these correlations, and long after the asymmetry in 
the evolution equations has become negligible. The novel 
droplet scaling state coarsens with I « 0.02£ in our re- 
duced physical units; this linear scaling is characteristic 
of the VH regime as in d = 3, but with a reduced prefac- 
tor. In the inert ial regime the system crosses over to a full 
symmetric morphology with a scaling law of I « 0.9 t 2 / 3 , 
similar to the d = 3 IH scaling state The ultimate 
fate of the IH regime in both d = 2 and d — 3 remains 
uncertain p9|,p[. 

Finally, we also performed some preliminary simula- 
tions of droplet morphologies for symmetric d = 3 fluids. 
As predicted in jnj, these morphologies always crossed 
over rapidly to a bicontinuous morphology, even when 
initialized in the VH regime: this contrasts strongly with 
our d = 2 results. Thus it remains possible that for 
d = 2 the absence of scaling of the bicontinuous mor- 
phology, and its restoration for a droplet state, is con- 
nected with the fact that the 50/50 bicontinuous state 
is (by symmetry) precisely at the percolation threshold. 



It would be interesting to see whether d = 3 systems 
close to percolation show similar anomalies, and whether 
a self-sustaining CIC droplet scaling is again recovered 
for compositions in that region. 

This work was funded in part by grant EPSRC (UK) 
GR/M56234. We thank I. Pagonabarraga for discussions. 



[1] A.J. Bray, Adv. Phys. 43, 357 (1994) 

[2] V.M. Kendon, J.C. Desplat, P. Bladon, M.E. Cates 
Phys. Rev. Lett. 83, 576 (1999); V.M. Kendon M.E 
Cates, J.-C. D esplat, I. Pagonabarraga, P. Bladon, cond- 



mat/0006026l J. Fluid Mech., in press. (The l(t) scaling 
plot was first used by S.I. Jury, P. Bladon, S. Krishna, 
M.E. Cates, Phys. Rev. E 59, R2535 (1999)) 
[3] A. Bhattacharya, S.D. Mahanti, and A. Chakrabarti, 

Phys. Rev. Lett. 80, 333 (1998) 
[4] Of course, these two morphologies have different mean 
density of the Euler invariant, and in that sense, as well 
as in the sense of percolation, have 'infinite range cor- 
relations'. Unlike correlations of the order parameter 0, 
however, these can evolve from a white-noise initial state, 
by a purely local early-time dynamics in finite time. 
[5] A.J. Wagner and J.M. Yeomans, Phys. Rev. Lett. 80, 

1429 (1998) 
[6] E.D. Siggia, Phys. Rev. A 20, 595 (1979) 
[7] H. Furukawa Phys. Rev. A 31, 1103 (1985) 
[8] The interfacial/inertial balance may be valid only at large 
length scales; since the interface area and kinetic energy 
are both decreasing, this cannot be true at all scales. See: 
V.M. Kendon Phys. Rev. E 61, R6071 (2000) 
[9] H. Furukawa, Phys. Rev. E 61, 1423 (2000) 
[10] H. Tanaka, Phys. Rev. Lett. 72, 1702 (1994) 
[11] V.S. Nicolayev, D. Beysens, and P. Guenoun, Phys. 
Rev. Lett. 76, 3144 (1996), V.S. Nicolayev, D.A. Beysens, 
Phys. Fluids 9, 3227 (1997) 
[12] R. Ahluwalia, Phys. Rev. E 59, 263 (1999) 
[13] A.J. Wagner, L. Giraud, and C.E. Scott, Comp. Phys. 

Comm. 129, 227 (2000) 
[14] Dynamic asymmetry can also arise through a <fi- 
dependent viscosity; but this continues to affect the equa- 
tions of motion themselves for all t. See K.E. Novik and 
PV. Coveney, Phys. Rev. E 61, 435 (2000) 
[15] E. Orlandini, M.R. Swift, and J.M. Yeomans, Euro. Phys. 

Lett. 32, 463 (1995) 
[16] M.R. Swift, E.Orlandini, W.R. Osborn, and J.M. Yeo- 
mans, Phys. Rev. E 54, 5041 (1996) 
[17] A.J. Wagner and J.M. Yeomans, Phys. Rev. E 59, 4366 
(1999) 

[18] For an animation of the phase-ordering please chec k 



http:/ /www. ph. ed.ac.uk/cmatter/links/awagner- drops 
[19] M. Grant and K.R. Elder, Phys. Rev. Lett. 82, 14 (1999) 
[20] F.J. Solis and M.O. de la Cruz, Phys. Rev. Lett. 84, 3350 
(2000) 

[21] U. Frisch, in Turbulence, Cambridge University Press, 
(1995). 



4 



